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BDS  Second-Order  NLTF 


I.  INTRODUCTION 


3 

Moat  circuits  in  a  typical  Cl  communication  system  are  nonlinear  to 
some  degree.  Examples  Include  preamplifiers,  mixers,  frequency-converters, 
channel  paths  containing  metal-to-metal-oxide  junctions,  and  in  particular 
power  amplifiers.  In  addition  to  these  inherent  nonlinearities,  there  may 
be  nonlinear! ties  deliberately  introduced  for  the  purpose  of  minimizing  the 
effects  of  the  inherent  ones.  These  circuits  usually  fall  into  a  class 
which  may  be  described  as  "mildly  nonlinear"  [»J]  circuits.  Since  these 
circuits  generally  have  memory,  a  simple  power  series  characterization  is 
usually  inadequate.  However,  a  Volterra  series  expansion  Cl]“[6],  which  is 
a  generalization  of  the  power  series,  provides  a  very  versatile 
characterization  of  a  nonlinear  circuit,  subsystem,  or  an  entire  system  [9] 
[14].  Furthermore,  the  Volterra  characterization  is  compact  for  mild 
nonlinearities  in  the  sense  that  a  truncated  Volterra  series  can  adequately 
describe  both  the  amplitude  and  memory  behavior  of  the  system. 

To  the  reader  familiar  with  Volterra  expansions,  the  basic  system 
entity  is  the  Volterra  kernel  h|^(T^  ,T2,  •••f  Tj^).  Its  Fourier  transform 
^k^^l ’^2’ *  * ’ ’ is  known  as  the  k-th  order  nonlinear  transfer  function 

3 

(NLTF)  [4], [8].  The  analyst  of  a  Cl  communication  system  (and  the  EMC 
engineer  responsible  for  the  design  and  implementation  of  the  nonlinear 
compensators  for  these  systems)  should  acquire  familiarity  with  techniques 
for  deriving  the  Volterra  NLTFs  from  the  circuit  or  its  equivalent 
description.  The  purpose  of  this  report  is  to  present,  in  a  practical  way, 
some  techniques  for  effective  representation  of  nonlinear  circuits/systeras 
by  Volterra  NLTFs. 


II.  BACKGROUND 


Numerous  alternative  representations  are  available  in  the  literature 
for  characterizing  and  analyzing  nonlinear  electronic  systems.  Of  these, 
the  Volterra  nonlinear  transfer  functions  (NLTFs)  description  [1],[2]  is 
particularly  attractive  since  it  lends  itself  to  convenient  frequency-domain 
interpretation.  As  such,  it  enables  straightforward  computations  of  such 
quantities  as  a)  linear  and  higher  order  nonlinear  responses  [2],  b) 
harmonic  distortion,  c)  intermodulation  distortion  [7],  and  d)  cross¬ 
modulation  distortion  [7].  Recent  research  has  shown  that  these  NLTFs  are 
also  well  suited  for  compensator  design  [15], [16]  to  minimize 
intermodulation  effects.  In  order  to  familiarize  the  reader  with  this 
analytical  and  design  technique,  this  study  briefly  introduces  the  Volterra 
expansion  and  then  uses  this  expansion  to  analyze  a  series  of  nonlinear 
phenomena . 

To  introduce  the  analytical  technique,  consider  the  input-output 
relationship 

y(t)-  T[x(t)]  (1) 

where  T  is  the  system  operator.  This  study  will  be  restricted  to 
relationships  which  are  time- invariant  and  only  "mildly  nonlinear."  For 
such  systems,  the  output  may  be  expressed  as, 

00 

y(t)=  I  y  (t)  (2) 

l<=1 

00  00  00 

=  I  /•  • '/hj^C  ,  •  •  •  ,  T|^)x(  t-T^  )  •  • 'XC  t-T|^)dT^  •  • ‘dli^ 
k=1  -“  -® 

where  y|^(t)=H|^[x(t)]  is  referred  to  as  the  k-th  order  response  and  is 
referred  to  as  the  k-th  order  system-operator.  These  various  notations  are 
consistent  so  that 

y^(t)»  H^[x(t)]  (3) 


/•••/  h,  (t.  .•••.T,  )x(t-T.  )'*'X(t-T,  )dT.  •••di, 


.-N  .  - 
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v** 


_l»  •  -  'Jl  •  ^  ^  - 


This  expansion  can  also  be  described  diagranunatically  as  appears  in 


Fig.  1. 


x(t) 


y(t) 


Fig.  1  Volterra  System  Representation 


This  expansion  of  y(t)  was  originally  described  by  Vito  Volterra  and 
later  named  the  Volterra  expansion  by  Wiener  [2]  who  applied  it  to  nonlinear 
noise  problems.  It  is  analogous  to  a  power  series  expansion.  As  with  a 
power  series  expansion,  this  "Volterra  expansion"  is  practically  useful  only 
if  the  series  converges  quickly  as  k  increases.  For  the  midly  nonlinear 
relationships  of  interest  in  this  report,  only  the  first  three  responses  h^ , 
h^,  and  h^  are  considered  significant. 

While  the  overall  nonlinear  relationship  of  equation  (1)  is 
nonhomogeneous,  equation  (3)  reveals  that  there  is  a  simple  relationship 
between  the  input  and  output  of  the  individual  k-th  order  responses  when  the 
input  is  scaled  by  a  constant  t.  Specifically, 

yk^(t)=  H^[ex(t)]  (iJ) 

=  e"  H^[x(t)] 

Thus  by  scaling  the  input,  the  factor  e  will  appear  as  a  multiplier  and 
can  be  used  to  identify  the  order  of  a  particular  response  [1].  This 
observation  will  be  useful  in  the  subsequent  sections  of  the  report. 


The  time-domain  integration  associated  with  these  expansions  are 
operationally  complex.  This  complexity  can  be  alleviated  by  use  of  the 
Fourier  or  Laplace  transformations.  In  the  image-space,  convolution  is 
isomorphic  to  multiplication.  To  demonstrate  this  fact,  and  to  determine 
the  proper  product  form,  a  multi-dimensional  response  y(|()  ^^1 '*'2* 
be  postulated 


00  00 

—  00  —00 


(5) 


where  is  referred  to  as  the  k-th  order  associated  response  [^4].  It  is 

apparent  that  this  associated  response  reduces  to  if  t^  =t2=  ••  •  =  t|^=t. 

But  the  associated  response  is  simple  to  Fourier  transform  to 

Then  Y,,(t)  is  simply  the  inverse  Fourier  transform  of  y/i.\ (i*.  •  *  "i  i",.) 

with  =  »  •  =  or 

y^(t)=  (6) 


•  •/{Hk(  Vi  ,  •  •  • ,  Vk)X(  '^l 


+  v.  )t 
k  dv. 


•dv, 


Finally,  the  Fourier  transform  of  Y|^(t)  becomes 

Y|^(f)=  (7) 

00  00 

•  ‘/{H  ( ,  •  •  • .  V  )X(  V  )  •  •  "Xi  V,  ) )  6(  f-(  V  +  •  •  •  ^v  ) )  dv  •  •  -  dv, 

-00  ^  '  <1  k  1  K  1  k 

This  equation  allows  easy  determination  of  the  k-th  order  response. 

i  2  ^  i  2  IT  F  t 

For  example,  suppose  that  an  input  x(t)  equal  to  e"^  1  +e'^  2  is  applied 

to  a  second  order  system  H2(f^if2)'  Then  equation  (7)  yields 

Y2(f)=  H^iF^ ,F^ )6(f-2F^  )  ^H^CF^ , F2)6( f-F^ -F^) 

+H2(F2,F^ )5(f-F^ -F^)  *H^(F^,F^)6{T-2F^) 


or  in  the  time-domain 


u  fc  I?  ^  .  j2’ir(F  +F  )t 


+  H2(F2,F^ 


In  general,  H^CF^.F^)  may  not  equal  H^CF^.F^).  But  often  it  is 
convenient  to  have  functions  which  have  this  "symmetry"  so  that 
H|^(f^  ,  •• -.f^)  equals  with  all  possible  permutations  of  the  independent 
variable.  This  can  be  guaranteed  by  defining  a  symmetrized  as 


—  H  (f 
k!  '  k^  1 


(8) 


where  the  tilde  indicates  that  the  function  is  unsyrameterized  and  the 
script-p  'P  denotes  the  summation  of  the  over  the  k~factorial 

permutations  of  the  independent  variables  [1]. 

In  the  above  example,  we  used  ,3^)  in  the  abstract  form.  A 

particular  realization  (although  not  the  most  general  one)  is  shown  in  Fig. 


y(t) 


Fig.  2  A  Simple  Second-Order  Volterra  System 


2.  Note  that  each  of  the  blocks  H  ,  and  H  is  linear.  For  this 

a  b  c 

structure,  it  can  be  shown  that 

H  (s, ,3  )=  H  (s  )  H  (s  )  H  (s  +s  ). 

212  a1  b2c12 

If  H  equals  H  ,  then  the  block  diagram  of  Fig.  2  can  be  more  concisely 

3  D 

depicted  as  in  Fig.  3- 


i**  JT,  .r - 
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x(t) 
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Fig.  3  Compact  Representation  of  a  Symmetric  Second-Order  System 


vv  .*  ^ 


III.  OVERVIEW  OF  EXAMPLES 


The  purpose  of  this  report  is  to  provide  techniques  for  deriving  the 
Volterra  NLTFs  from  other  available  descriptions.  To  this  end,  the 
subsequent  sections  will  introduce  the  method  of  analysis  through  a 
collection  of  nonlinear  electronic  device  and  circuit  examples.  These 
"available  descriptions"  constitute  a  mix  of  differential  equations,  device 
representation,  circuits,  and  system  representation  by  block  diagrams. 

The  first  two  examples  consider  nonlinear  differential  equations.  The 
first  example  examines  a  nonlinear  second-order  differential  equation;  the 
coefficient  of  the  first  derivative  is  a  variable  and  dependent  upon  the 
output.  The  second  example  is  a  generalization  of  the  first  equation;  it  is 
an  n-th  order  differential  equation  where  all  of  the  coefficients  are 
variable  and  dependent  upon  the  output. 

The  third  example  analyzes  a  nonlinear  device  wherein  the  current  can 
be  expressed  as  a  power  series  expansion  of  the  voltage;  a  circuit  including 
this  device  is  analyzed  using  nodal  analysis.  A  forward  biased  diode  has  a 
nonlinear  current-voltage  relationship  as  described;  but  the  diode  is 
further  complicated  by  the  voltage-dependent  capacitance  associated  with  its 
junction.  Therefore,  example  four  analyzes  the  diode  as  a  device  outside  a 
circuit.  This  example  is  particularly  important  because  of  the  wide  spread 
use  of  diodes,  and  because  the  nonlinear  characteristics  of  diodes  can  be 
employed  in  circuits  designed  to  act  as  "nonlinear  compensators." 

The  fifth  example  is  the  complementary  dual  to  example  three.  Here, 
the  nonlinear  devices  have  a  current-voltage  relationship  such  that  the 
voltage  can  be  expressed  as  a  power  series  expansion  of  the  current.  Also, 
the  analysis  employs  loop  equations  rather  than  node  equations. 

The  sixth  example  employs  a  more  complicated  nonlinear  device,  the 
transistor.  Here  the  transistor  current  is  modeled  as  the  "product-power" 
series  expansion  of  two  voltages. 

The  seventh  example  analyzes  a  cascade  of  two  nonlinear  systems.  This 
example  is  particularly  important  because  this  is  a  configuration  which  can 


be  employed  to  eliminate  the  nonlinear  output  of  the  first  system.  The 
example  explicitly  describes  how  to  develop  a  nonlinear  cascade  compensator 
to  eliminate  the  nonlinear  response  of  the  first  system.  The  eighth,  and 
final,  example  uses  the  nonlinear  characteristics  of  a  diode  to  implement  a 
nonlinear  cascade  compensator. 


Volterra  nonlinear  transfer  functions  can  be  employed  to  characterize, 
in  the  frequency-domain,  certain  classes  of  nonlinear  differential 
equations.  As  a  specific  example,  consider  the  equation 

— 2  ^  dt  f(y)}  x(t)  (9) 

dt 

If  f(y)  can  be  expanded  in  a  power  aeries  which  converges  quickly,  then  y 
may  be  expanded  in  a  Volterra  expansion  which  also  converges  quickly. 


Solution. 


It  was  stated  that  f(y)  can  be  expanded  as  a  power  series;  i.e., 


f(y)-  I  a_y" 

n»0 


Now  form  a  Volterra  series  expansion  for  y  as 


y=  H[x]=  I  H  [X]-  ly  (11) 

k=1  k-1 

Substituting  these  expansions  into  the  differential  equation,  one  obtains 

-2  ^  ^  ^k^"  ^  ^  ^  ^  ^  ^  ^n  ^  ^  ^  ^k^  “ 

dt'^  k=1  j=1  ^  n=0  k=1  k=1 

For  this  example,  the  first  three  Volterra  NLTFs  for  y^  ,  y^,  and  y^  will  be 
formed.  This  is  accomplished  by  explicitly  writing  the  individual  terms  of 
the  various  orders  for  the  above  differential  equation.  To  keep  track  of 
the  order,  scale  the  forcing  function  x  to  ex. 


e-Scaling^ 


^e  “  ^  Hj^Cex]  =  I  e  H|^[x] 
k»1  k=1 


.*•  ,  “  . *■  "  >  »*-  «"*• 
’*•*•*“•*  •’  •  ii*  ^ 


The  differential  equation  then  becomes 


— P  {ey, Yp+e 

dt""  ^ 

"  ^  {aQ>ey^^a^e2y^2^e3jaQy3-2a^y^y2-^a2yi^}--' 

2  3 

+b{£y^+e  y2‘*’e  y^"^***}  * 

Now,  like  powers  of  e  can  be  collected  and  equated. 


S  !'l  *^0  ^  'l  •‘■i'l  ■  * 

dt 


Or  more  concisely, 


G[y^]=  X, 


d^  d 

where  G  is  the  linear  differential  operator  {— ^  +3^  ^  y^  +b} 

dt 


y^-  G  '[x]  =  H^[x]. 


This  equation  can  be  Laplace  transformed  to 


-2 - : 

s  +aQS+b 


It  therefore  follows  that  is  G  or 


H^(s)=  - 

s'^+aQS+b 


d  d  w  «r  n  d  2 
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dt 
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So  that 


■  3t  »,[*])] 


Now  the  Laplace  domain  formulation  of  can  be  performed  [U  ].  The 
determination  of  this  NLTF  is  particularly  simple  if  is  diagrammed  as 
described  in  Section  II. 

”^1^1 


^2  h  **0  Ht  '3  "’'s"  '  it  Wi^> 


and  therefore 

.,-1  rd 


Here  again,  H,  may  be  more  apparent  from  the  block  diagram  of  Fig.  i(. 


Fig.  4  Block  Diagram  of  the  Third-Order  NLTF 
From  the  diagram,  H_  can  be  written  as 


where  the  tilde  has  been  included  over  to  indicate  that  this  is  an 

unsymmetrized  transfer  function.  The  symmetrized  is  obtained  from  as 

described  in  Section  II  as  H^-  ^  H^Cs^.s^.s^)  [1],  And  therefore,  the 

symmetrized  is 


-H^  (s^  +S2  +  32){3^ 


2a 


+  {H^  (s^  )H2(s2,s  )  +H^  (s2)H2(s^  ,s  )  +H^  (s  )H2(s^  .s^)} } 


Here  it  is  noted  that  the  y^^'s  a>"e  recursive  in  the  sense  that  y^^  is  only  a 


function  of  y^  ,  y^,  •••,  y|^_^  •  Thus  the  truncation  of  the  original  equation 


to  the  third  order  terms  leads  to  no  error  in  the  development  of  the  first 
three  NLTF's.  Alto,  the  process  could  be  continued  to  obtain  as  many  of  the 
NLTF's  as  required.  Indeed,  there  Is  no  mathematical  reason  that  all  the 
y’s  could  not  be  formed.  But  in  most  applications,  this  Volterra  expansion 
is  only  practical  for  analysis  and  design  if  only  the  first  few  terms  are 
significant. 


Example  2;  General  Differential  Equation. 


The  first  example  was  a  differential  equation  which  was  "nonlinear  in 
the  second  term."  This  is  a  special  case  of  the  nonlinear  differential 
equation 


d 

dt 


.n-1 


n  {  I  a  y  }  +  —  {  I  a  y  ^  a  y  }  = 


(19) 


This  is  a  rather  general  nonlinear  differential  equation  since  it  can 
represent  a  n-th  order  differential  equation  where  each  coefficient  of  the 
differential  terms  are  dependent  on  the  output,  but  can  be  expressed  as 
power  series  expansions.  The  approach  to  solving  this  equation  is  identical 
to  the  previous  example. 


Solution^ 

Expand  y  as  a  Volterra  series 
00  00 

y=  H[x]=  I  H  [x]  -  I  y^. 
k=1  k-1 

As  before  this  expansion  is  substituted  into  the  differential  equation  and 
the  input  or  forcing  function  x  is  scaled  by  e  to  keep  track  of  the  order. 


e-Sc§Hng^ 


d 

dt 


n  <» 


n  t  I  a  {  I  {  [  a  {  I  ♦••• 

"  11=1  k=i  dt  ’  £=1  k=1 


n-1  ® 


k 


I  %  J  I  cx 

£=1  k-1 


This  equation  can  be  rewritten  as 


Now,  by  defining  the  linear  operator  as 


n  .i 

V  d 


4“  ^  -i  ^  H 

i-1  dt  ’ 


the  differential  equation  (20)  can  be  rewritten  as 


00  OO  00  )<+•.. >)< 

I  1..  I  ...  I  e'  '■ 


■t  "r’  “t" 


Now  as  before,  this  expression  will  be  evaluated  for  various  powers  of  e, 
and  like  powers  will  be  collected. 


^CYi]-  X, 


so  that  is  ;  or  in  the  complex- domain 


H^(3)  = 


n  n-1  ^ 

a  , s  +a  ,  , 3  + • • •  *3^  , 

n,1  n-1,1  0,1 


L^Cy^]-  -L^Ly,^] 


so  that 


y^-  -hJL2[H^[x]H^Cx]]] 


The  Laplace  domain  formulation  of  produces 
,3^)”  -H^  (3^  )H^  (s2)L2(s^ +32)H^ 


•*.  A  •'  •*.  •* 


?  7  J*  ^ 
V  S.  %  ' 

^ 

■  »  •  «_•  •  • 

Xl>CvW 


■  -•  *>.*  ' 


m 


•<*  <■ ' 


'  V  .* 

■  ■  V 


vN-’-v 
V  :y\' ' 


L^H^-  -2L2[H^H2] 


then  the  Laplace  transform  is 

.s^.s^)-  -2L2(3^+32+32)H^  (3^+32*32)H^  (s^  )H2(S2+s^)  (24) 

-  L2(3^+32^3^)H^  (s^+32+32)H^  (3^  )H^  (S2)H^  (3^) 

where  the  tilde  has  been  included  as  a  reminder  that  this  is  an 
unsymmetrized  transfer  function.  The  symmetrlzation  of  this  transfer 
function  is  straight  forward  as  described  in  Section  1  and  as  applied  in 
example  1 . 

4 

£  ; 

LiCH4]=.  -L2[H2H2+2H^H^]  ] 

So  that  Hjj  can  be  formed  in  the  complex-domain  directly  or  from  a  diagram  as 
H^(3^  ,32,32. Sj^)-  (25) 

-L2(s^+32+32  +  3^j){H2(3^  ,S2)H2(32,3jj)+2H^  (s^  )H2(s2,S2,3^)} 
-3L2(3^+S2>32  +  3j^)H^  (s^  )H^  (s2)H2(32,3^) 

“S  ^  ®2*  ®4  ^  ^  ®2  ^  ®3  ^  ®4  ^ 

As  stated  earlier,  these  transfer  functions  generally  are  recursive  so 
that  the  complexity  of  the  function  grows  with  the  order.  This  is  apparent 
in  this  example.  It  is  for  this  reason  that  this  analysis  is  usually 
applied  only  to  slightly  nonlinear  systems  where  it  is  likely  that  the 
Volterra  series  expansion  will  converge  rapidly. 


15 


The  objective  of  this  analysis  is  to  find  the  form  of  the  Volterra  NLTFs  H 
which  relate  the  node  voltage  to  the  input  voltage;  i.e., 


y-  I  H^[v]-  I  y 

k-1  k-1 


Solution. 


If  the  input  v  is  scaled  to  ev,  then  a  new  node  voltage  y^  results 

ao  eo  00 

y^-  I  I  H|^[v]-  I  y^^ 


Substituting  this  into  the  node  equation, 


(28) 


^  d  p  d  r  k 

e  C  —  V-  C  ) 
dt  dt,  ‘* 


0»  00 


k  in 


1  'k  *  'C  b,  ‘"'k  *  Vn^wV'^'ki 

1  •'k-1  n-1  k-1 


00  00 


g[  I  e'^y^]  I  SrI  I 

k-1  n-2  k-1 


00  00  00  00  00  )(  -f  •  •  •  4')< 

G[  I  *  I  \  I  I  ••*  Z  e  ^  ^  "  ^k  ^k  "‘^k 

k-1  n-2  k-1  k.-l  k  -1  1  ^2  n 

1  2  n 


where  G  is  the  linear  operator  with  the  Laplace  transform  Cs+a  +  — ;  i.e., 

1  LiS 

G(s)  is  the  linear  admittance. 


Now  as  stated  earlier,  we  are  only  interested  in  the  first  three 

1  2 

NLTFs,  and  therefore  only  terms  that  contribute  to  e  ,  e  ,  and  e  need  be 
explicitly  expressed  in  the  equations  (28).  So  equation  (28)  can  be 
expressed  as 

eC  ^  (29) 

e^G[y,]  +e^lGCy2]  ^e^CGCy^] 

where  HOT  are  the  higher  order  terms  in  e.  Notice  that,  this  finite 
expansion  of  equation  (29)  is  more  readily  obtained  from  the  double  sum  form 
of  equation  (28)  rather  than  the  last  form.  While  the  last  "infinite  sums" 
form  of  equation  (28)  is  consistent  with  the  general  theory  of  Vol terra 
functions,  for  this  example,  and  for  all  the  future  examples,  this  infinite 
sums  form  is  operational  less  efficient  then  the  double  sum  form. 

Therefore,  in  the  future,  this  more  complicated  form  will  not  be  employed. 

Using  equation  (29),  the  terms  associated  with  the  powers  of  e  can  be 
collected. 

1 

E  ; 

C^v.  GCy,]  (30) 

This  equation  can  be  Laplace  transformed;  and  then  by  defining  a  linear 
impedance  Z(3)  as  G  ^(s),  y  becomes 


(3)-  Z(3)  Cs  v(s)  (31 ) 

and  therefore 

H^(s)-  Cs  Z(3)-  Cs  - -  (32) 

LCs'^  +La^  3  +1 


2 

e  : 

G[y2>  (33) 

The  Laplace  transform  of  can  then  be  formed  from  the  above  equation,  or  a 
block  diagram  of  the  above  equation,  as 

{H^  (s^H^  (s^)}  (3^) 

J. 

e _ 

GCy^]^ 

And  therefore, 

H^Cs^  ,32,3^)=  (36) 

-ZCs^+s^^s^) |2a2H^ (3^ )H2(S2,s^)  ^a^H^ (s^  )H^ (32)H^ (s^) } 

where  the  tilde  has  been  included  over  as  a  reminder  that  this  expression 

Is  unsymmetrized  .  The  symmetrized  is  obtained  from  as  described  in 

Section  1  as  H2=  -jj-  ■  H(s^,32,3^).  The  symmetrized  e  .pression  is 

therefore, 

H^(s^  ,82,3^)=  (37) 

-Z(3^+32  +  S2)  |a^H^  (s^  )H^  ^®3^ 

+(32/3) {H^ (s^  )H2(s2,s^)+H^ (32)H2(s^  , )+H^ ( )H2( s^  .S2)}} 


These  three  NLTFs  can  be  depicted  graphically  as  appears  in  Fig.  6. 


The  previous  example  Included  nonlinear  devices  such  that  the  current 
through  them  can  be  expanded  as  a  power  aeries  of  the  voltage  across  them. 
One  example  of  this  type  of  nonlinear  device  Is  a  forward  biased  diode.  But 
the  diode  Is  more  complicated  since  there  are  capacitances  and  resistances 
also  associated  with  It.  The  capacitance  la  a  function  of  thickness  of 
the  junction  regions,  and  this  thickness  Is  Itself  a  function  of  the  voltage 
across  the  diode;  so  that  the  capacitance  Is  not  a  constant  but  depends  upon 
the  junction  voltage.  Since  the  diode  la  such  a  common  nonlinear  electronic 
device,  It  will  be  analyzed  In  this  section  as  an  Isolated  device. 

As  stated  above,  a  power  series  representation  will  be  used  to 
describe  the  voltage-current  relationship  of  the  diode.  The  order  of  this 
power  series  expansion  generally  depends  upon  the  quiescent  point  and 
operating  range.  In  many  practical  applications,  the  operating  conditions 
are  selected  so  that  only  the  dc,  linear,  and  quadratic  terms  are 
significant.  But  for  the  present  analysis,  the  cubic  term  will  also  be 
included. 


+ 


y 


(a) 


(b) 

Fig.  7  Diode  Models 


(c) 


-  20  - 


r- 


•  . 


Fig.  7(b)  depicts  a  standard  model  for  a  diode.  Typically,  R„  is  of 

D 

the  order  of  a  few  tenths  of  an  ohm,  and  R,  is  several  hundred  kilo-ohms. 

L» 

For  the  applications  of  interest  to  this  study,  these  resistances  are 
inconsequential.  Therefore  the  simpler  model  of  Fig.  7(c)  will  be  employed 

The  objective  of  this  analysis  is  to  determine  the  first  three 
Volterra  NLTF's  and  P^  which  relate  the  voltage  to  the  current. 


Solution. 


If  the  diode  of  Fig.  7(c)  is  driven  by  a  current  source  x(t),  a 
voltage  y(t)  will  result.  Furthermore,  the  nonlinear  current  source  and 
capacitance  of  the  diode  model  will  have  the  forms 


‘d- 

n-1 


where  the  voltage  y  can  be  expanded  as  a  Volterra  series 


y-  I  PJx]  -  I  y, 


The  node  equation  for  the  diode  and  external  current  source  is 


dt  ° 


d  V  n+1  V  n 

-TT  1  c  y  ■‘■/ay  -x 

dt  n-'  n'' 

n»0  n-1 


1  ^  -* 


e-Scaling^ 


Now  the  driving  current  source  x  can  be  scaled  by  e  to  obtain 


--if 


(39) 


Vvi  3t  *vti J ° 

n»l  K*l 


^  ^  ^a2}[e^y,^  +2e3y^y2] 


+{C2  ^  +a2}[e^y^3  -ex  +HOT 

where  HOT  denotes  the  higher  order  terms  in  e.  The  powers  of  e  can  now  be 
collected  and  equated. 

1 

e  : 


Here  it  is  convenient  to  introduce  the  notation  L  (s)«  c  , s  +a  .  Then  the 

n  n-1  n 

Laplace  transform  of  the  above  equation  becomes 
(s)  y^ (s)=  x(s) 

From  this  equation,  the  Laplace  transform  of  is  found  to  be 


(40) 


2 
e  : 


Fig.  8 


Second-Order  NLTF  of  a  Diode 


The  form  of  the  NLTF  is  apparent  from  Fig.  8  as 


(41) 


1 

\ 

I 


e3: 


”^°1  dt  ^2^  ~^°2  dt  ^ 

-  -h^l2Y^y^'\  -L^Cy^^] 

The  multidimensional  Laplace  transform  of  P^  can  then  be  formed  as 


P^(s^  .s^.s^)-  -P^(s^+S2+32){2L2(s^+S2+32)P^(s^  )P2(32,3^) 

+  L2(s^+S2+32)P^  (3^  )P^  (32)P^  (s^) } 

where  the  tilde  has  been  included  over  P^  as  a  reminder  that  this  is  a 
unsymmetrized  operator  as  discuased  in  Section  1. 


(42) 


Thia  then  completes  the  analysis  of  the  diode. 


VI.  NLTFS  OF  MULTI-LOOP  AND  DEPENDENT  SOURCE  CIRCUITS 


Example  5-  Nonlinear  Multi-loop  Circuits, 


Fig.  9  Multi-Loop  Circuit  with  Nonlinear  Devices 


The  previous  two  circuit  examples  dealt  with  nonlinear  devices  such 
that  the  voltage  across  those  devices  could  be  described  as  a  power  series 
expansion  of  the  current  through  them.  For  this  example,  the  opposite  is 
true;  the  nonlinear  devices  are  such  that  the  current  through  them  can  be 
described  as  a  power  series  expansion  of  the  voltage  across  them. 

Fig.  9  depicts  a  network  with  two  static  nonlinear  devices  which  are 
represented  as  D^  and  D^.  Here  the  R's  represent  "linear  impedance 
operators"  and  actually  reflect  the  presence  of  any  linear  "sub-network." 

The  loop  equations  for  Fig.  9  can  be  written 

v=  ^  [x]  - 

0=  -R^^Ex]  +  ~  ""^2 

where  *^1 1  “**10^**1 2’  **22°**20^**1 2’  voltages  u^  and  u^  are 

nonlinearly  related  to  the  currents  through  the  devices  D^  and  D^ 
respectively. 

The  general  objective  of  this  analysis  is  to  relate  the  currents  x(t) 
and  y(t)  to  the  input  voltage  v(t).  But  more  specifically,  the  objective  is 


-  24  - 


to  find  the  form  of  Volterra  transfer  functions  P,  and  Q,  which  relate  the 

k  k 

T 

k-th  nonlinear  current  [X|^  to  the  Input  v(t);  viz., 


c(t)-  I  X  (t)-  I  P  [v(t)] 
k»1  k-1 


y(t)«  I  y^(t)-  I  Q  [v(t)] 
k-1  ^  k-1 


For  efficient  notation,  define  as  a  column  vector 


-k°  Q, 


Solution^ 

We  begin  by  requiring  that  the  voltages  and  can  be  adequately 
modeled  as  a  power  series  expansion  where  only  the  first  three  terms  are 
significant;  i.e.. 


Ui=  I  b  (x-y)"-  I  b  c" 
n-1  n-1 


n-1 

where  the  current  through  the  common  branch  is  C“X-y. 


e-Scalingi 

Now  replace  v(t)  with  a  scaled  voltage  ev(t)  and  thus  generate  new 
currents 


X  =  I  "  I  E  ^  ^  \ 

^  k-1  k-1  k-1 


y  “  I  Q^l^ev]  -  [  e  Q  [v]  -  Icy. 

^  k-1  k-1  k-1 

Substituting  this  back  into  the  loop  equations: 


*,  i*.  *.  •.  *. 


cv.  I  cS  -R,2[yx]l  •!  »„(J  'V" 

k=1  n=1  k»1 


0-  I  "«22‘^yk^f  I  -I  a  {  I  e\}” 

k»1  n-1  k-1  n-1  k-1 

As  stated  earlier,  we  are  only  interested  in  the  first  three  NLTFs,  and 

12  ^ 

therefore  only  terms  that  contribute  to  e  ,  e  ,  and  e  need  be  explicitly 
expressed  in  the  equations  (^6).  So  the  first  equation  of  (46)  can  be 
expressed  as 

ev.  j  cN  -i  b„(J  x\)"  .HOT 

k=1  n»1  k»1 

where  HOT  represent  the  higher  order  terms  of  e  as  k  runs  from  four  to 
infinity,  and  also  as  n  runs  from  4  to  infinity.  Expanding  the  summations 
and  regrouping  the  terms  results  in 


j*  t** '  •  *  • 


t*  ^ 

--v'Vv'.vi'.CS 


GV=  e’  {R^^[x^]  +b^ 


•^e^lR^^Cx^]  -R^2^y3^  *^1^3  *^^2^1  ^2 


where  HOT  now  includes  the  previous  higher  order  terms  and  also  the  higher 
order  terms  formed  from  the  square  and  cube  processes.  It  is  emphasized 
that  the  expressed  terms  are  exact.  The  above  manipulations  have  merely 
ignored  terms  of  e  which  have  power  greater  than  three. 


.vT'.y. 


The  second  equation  of  (46)  can  be  manipulated  to  result  in 

-e^I-R^^Ex^]  -^,^2  ‘^2^1^  "^1^2  "Vl^^ 

+  e^l-R2^[x,]  ”*^3^1  "’®1^3*^®2^1^2'"^3^1 


Now  the  terms  associated  with  the  powers  of  e  can  be  collected  and  equated. 


0-  *«22Cy,]  -!>,«,  ‘^y, 

or  more  concisely, 


V 

=  it- 

(R^^+b^) 

*1 

0 

-(R2^+b^)  (**22*^1  “^^1^ 

,^1 

[Z] 


where  [Z]  represents  the  "linear  network  impedance." 


(1*9 


Taking  the  Laplace  transform  of  equation  (49)  yields 
\j/(s)=  [Z(s)]  _i^  (s) 


The  "admittance  matrix"  CG(s)]  can  now  be  formed  as  the  inverse  of  the 
matrix  [Z(s)]; 


[G(3)]= 


[Z(s)] 


1 

A(s) 


(R2,(s)+b^ ) 


(R^,(s)+b^ ) 


(50 


where  A(s)  is  the  determinant  of  the  matrix  [2(s)].  Then  it  follows  that 
(S)»  [G(3)]  Jj;(s) 


1 

A(3) 


v(s) 


But  since  operating  on  v(t)  is  defined  as  i^ ,  it  follows  that 

(R22(s)'^b,  *^1  )“ 

(R,2(3)-bi)  _ 


H^(s)= 


P^(s) 

1 

Q,  (3) 

”  A(s) 

(51 


The  relationship  between  _i^  and  v  can  be  depicted  graphically  as  appears  in 
Fig.  10. 


H^(3) 


(t) 

y^Ct) 


Fig.  10  First-Order  NLTF 


0-  *^2^1 

0-  -S^iCXj]  '"a/yp  -“,52  -“2^1^  •“l'2  *Vl^ 

or  more  concisely, 


*2  ^2^1 

0=  [Z]  /  >  2,2 

^2  2^1  2^1 


V  •»'  '.*  •/  • 

■'•*.%  A  ■.'•V 
S'.VC.N-. 
-.'.vV-'A'. 


i.-  CG]  ,%'2 


"2^1  2''1 


This  relationship  can  be  depicted  graphically  as  appears  in  Fig.  11 


---'1^  v>^-.>cV4  ~  ^ 

Fig.  11  Second-Order  NLTF 

But  to  solve  for  the  form  of  and  Q^,  equation  (52)  should  be  Laplace 
transformed.  Since  these  equations  Involve  a  squaring  process,  the 
multidimensional  Laplace  transforms  are  employed.  For  convenience,  define 

M  =  P  -Q 

k  k  k 

Then  the  transfer  function  from  v  to  <(i  in  Fig.  11  is  -b^M^  ( s^  )M^  (s^) ;  the 
transfer  function  from  v  to  0  is  b2Mj (s^  )M^ (32)-a2Q^  (s^  )Q^ (s^).  And 
therefore,  the  transfer  function  from  input  to  output  is 


^1  '' 
— "V  2 


CG(s)J 


•  V 


t  ■'M.  *-• 


• .  • 


H2(3,.32)=  [GO^^-S^)] 


-b^M^ (3^ )M^(32) 


_b2M^(3^  )K^(32)  -a^Q^  (3,  )Q/32)_ 


0=  R^^[x3]  -R^2[y3]  *Vl' 


0-  -R„Cx_]  ^R^^Cy,]  -2b2C,C2  -b3C/  ^a^y3  ^  ^a^y^y. 


'21^*3"  22^-'S 

or  more  concisely, 


^3 

1 

0“  [Z] 

+ 

■i  3 

Jl 

l::S'‘2'^t«2  *^Vi>'2  *Vi- 

and  therefore, 


13=  [C] 


JbjtiSj  »b^5,^  -^W2  'Vi- 


This  relationship  can  be  depicted  as  appears  in  Fig.  12. 


Fig.  12  Third-Order  NLTF 


The  transform  function  from  v  to  and  Is  apparent  from  equation  (54)  as 
H^Cs^.s^.Sj)  (55) 

-2b2M^  (s^)M2(S2,s^)  -b^M^(spM^(S2)M^(s^) 

-  CG(3^+S2+s^)] 

2b2M^ (s^ )M2(s2,S2)  +b^M^ (s^ )M^ (32)M^ (S^) 

-2a2Q^ (s^ )Q2(32,32)  -a^Q^ (s^ )Q^ 

where  the  tilde  has  been  included  over  as  a  reminder  that  this  expression 
is  unsyimnetrized.  Again,  the  symmetrized  is  obtained  as  described  in 

Section  1  as  ‘  stated  objective  of  this  example 

was  to  determine  the  first  three  P's  and  Q’s  and  that  has  been  accomplished. 


Example  6;  Dependent  Sources  (Transistor  Model). 


This  example  analyzes  a  circuit  with  a  nonlinear  transistor  model. 
Here  the  transistor  current  is  modeled  as  a  function  of  two  voltages  u  and 
w;  viz. , 


T  V  V  n  n 

1=2.  2.  g  u  w 

^  nm 

n=0  m»0 


where  gQQ=0*  The  circuit  of  interest  appears  below  in  Fig.  13* 


w  + 


Fig.  13  Transistor  Model 


Generally,  the  objective  of  this  example  is  to  determine  the  output 
voltage  y  as  a  function  of  the  input  voltage  x.  More  specifically,  the 
objective  is  to  determine  the  first  three  Volterra  NLTF's. 


1  /K 
-1/R, 

V 

1 

X  ^ 

-1 

L_ 

I  I 

n»0  m-0 
n+m>1 


n  m 
g  u  w 
“ran 


where  the  linear  part  of  the  transistor  model  has  been  moved  to  the  left 
side  of  the  equation. 


G-Scaling^ 

The  voltages  u  and  w  can  be  expanded  as 


u=  P[£x]=  I  P|,[ex]  =  5;  e  u 


(58) 


k=1 


k=1 


w=  Q[ex]=  I  Q;^[ex]  =  I  w 
k»1  k*1 


where  the  scaling  constant  e  has  been  included  to  keep  track  of  the  order  of 
the  Volterra  operators.  These  expansions  for  u  and  w  can  now  be  substituted 
into  the  node  vector  equation;  viz., 
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1  /R 
-1/R, 

1 

X 
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L_ 

V  V  /V  ><  k  (59) 

n=0  m=0 


n+ni>1 


where 


Now  like  powers  of  e  can  be  collected  and  equated 
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1/R^ 

-1/R, 


or 
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-1 


1  /R^ 
-1/R, 


X  =  X 


where 
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«  d  1 
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1  dt  R^  *^01 
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^01 
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2  dt  ^R^  R^  ‘^10_ 


1/R^ 

-1/R, 


and  A  is  the  determinant  of  [L]. 
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(61  ) 


(62) 


and  therefore  H.  can  be  described  in  "operator  notation"  as 


«2“  ^802^1  ^ 


(63) 


(830^^  "82iU^S"8i2'^i“i^  "803“i^ 
+2g2oU^U2  *811  ('^i”2‘^'^2”i  ^  '^^8o2'^l”2^ 


Rj^-1  ^830^1  ^821 P,  ^i'*'8i2^iS  *8q3Qi^ 


■"^820^1  **2  '"81 1  (Pi  Q2‘'^2S  ^  "^^802^  *^2^ 


where  H^  is  expressed  in  operator  form. 

As  in  all  the  previous  examples,  is  recursive  in  the  sense  that  it 
is  a  function  of  only  previous  H*3. 

Now  the  operators  can  be  Laplace  transformed  to 

"Cs>l/R 

Hi(3)=T7^  (66) 

^  +11/R,+1/R2}_ 

-2^®1  ^®1*®2^^82oPi  (®i  )Pi  (®2^'^8i  1  P( ^i  )Qi  (®2^^8o2S  (®i  (®2^^' 

H3  ( ,  s^ ,  83 )=  ( s^ +  S2'*'33  )  (  g^Q?^  ( 3^  )P^  ( s^ )P^  ( 3^ ) 

■*■821  Pi  (^i  )Pi  (®2^®1 

^803^1  (®1  ^®3^  *^82oPi  (®i  )P2(®2’®3^ 

+g^^  {P^  (3^  )Q2(32,S3)'*-P2(3,  .32^'^1  ^®3^ 


*^802^  (®t 


where 


A(s)=  C^C23^  +{C2(gQ^+1/R^^)+C^  (1/R^+l/R2-giQ)}3 

^{(8o,-1/Rl)(1/Ri^1/R2)  -8io/Rl} 


_  3ii  _ 
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VII.  NLTFs  OF  CASCADED  SUBSYSTEMS 


Example  Nonlinear  Cascade. 


In  this  example,  we  wish  to  design  a  cascade  filter  or  compensator  P 
which  will  eliminate  some  of  the  nonlinearities  of  a  system  H,  as  depicted 
in  Fig.  li*. 
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Fig.  14  Cascade  of  Two  Nonlinear  Subsystems 

Here  it  is  assumed  that  the  form  of  the  nonlinear  system  has  already  been 
determined  as 

00  00 

W=  H[x]=  I  H  [x]-  I  w  (61 

n=1  n-1 

Indeed,  after  the  form  of  the  compensator  has  been  determined,  the  H's  for 
an  earlier  example  (viz.,  the  single  nonlinearity  example)  will  be  used  to 
Implement  part  of  a  specific  compensator. 

Solutlon_^ 

Form  the  output  y  as 
00 

y=  P[w]=  I  P,^[w] 
k^l 

00  00 

k=1  n=1 

Then  replace  x  with  cx  to  form  a  new  output; 


(69) 


J,  ’’kt  I 

k-1  n«1 


00  00 


■  I  pj  !  t\cx]) 

k-t  “  n-1 

This  relationship  can  be  reexpressed  in  terms  of  the  "x-y  transfer  function" 
Q  as 

J, 

Now  if  powers  of  e  are  collected,  the  Q's  can  be  determined  using  "operator 
notation"  as 


Pi”! 


(70) 


Q  =  PH  *P  H 
^2  12  2  1 


+2P2lH^.H2}  ^P^H^ 


P^H^  +2P2{H^.H2}  ♦P2H2  +3P3{H^,H^.H2} 

Q^*  P^H^  +2P2{H^,H^}  +2P2{H2,H2} 

+3P3{H,  ,H2.H2}  +4P^j{H^,H,,H,,H2}  ^P^H^ 


These  equations  can  be  expressed  in  the  complex-domain  as 
Q^(s)=  P^(s)H^(s) 

02(81  ,S2)=  P^  (s^+32)H2(s^  ,S2)+H^(s^  )H^(s2)P2(3^  .s^) 


P/Si.Sj.SjjHjCSi.Sj.Sj)  .  J  CH,(s,)H2(S2.33)P2(8,,82.Sj) 


+H^  (S2)H2(s^  ,s^)P2(s2,s^+s^)+H^  (s^)H2(s^  ,S2)P2(s3,s^ +82)] 


The  stated  objective  of  this  exercise  is  to  reduce  the  nonlinearities 
without  distorting  the  linear  characteristics  of  H.  Therefore,  the  form  of 
should  be  ;  and  therefore, 

Pl=  I  (71) 


where  I  is  the  identity  operator.  Furthermore,  and  can  be  made  zero 
by  choosing  and  as 

P2=  ^■^2) 

P^=  -(H^  +2P2H^H2)H^'’ 

Now  it  was  assumed  that  only  the  first  three  H's  were  significant.  So 
choose  the  remaining  P's  as  the  null  operators.  Then  the  remaining  Q's  are 
only  functions  of  the  insignificant  H's  and  "cross-products"  of  the  first 
three  H's.  Thus  by  including  only  the  first  three  P's  in  the  compensator, 
the  first  two  "high-order"  nonlinearities  are  eliminated,  and  the  remaining 
nonlinear  terms  are  small. 


In  order  to  be  more  specific  concerning  the  implementation  of  a 
compensator,  let  us  implement  the  P^  part  of  the  compensator  using  the  H's 
of  an  earlier  model.  Specifically,  the  H's  of  the  "single  nonlinearity 
example"  is  employed. 


H  (3)=  Cs  - - 

LCs'^  +Lg^s  +1 

H^CSi.s^)"  -Z(S^+S2)  ^2  {H^(s^)H^(s2)} 
Then  becomes 
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^  LC(s  +s  +Lg  (s  +s  )  +1 


•  A  k.*» 


W  H  '■ 
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The  algebraic  expression  for  P3(s^ . 82.82)  Is  determined  by  evaluating  the 
second  expression  of  equation  (72).  While  this  evaluation  is  mathematically 
straight  forward,  the  resulting  expression  is  involved.  Therefore,  the 
specific  form  of  will  not  be  determined.  This  then  completes  this 
example. 
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VIII.  NLTFs  OF  A  PROPOSED  COMPENSATION  NETWORK 


Example  8:  Diode  Compensator. 


The  previous  example  dealt  with  the  general  design  of  a  cascade 
compensator  to  reduce  nonlinearities.  This  example  analyzes  a  specific 
circuit  (viz.,  a  Balanced  Diode  Squarer  (BDS)  as  depicted  in  Fig.  15)  which 
can  be  employed  as  part  of  such  a  compensator.  Specifically,  it  can  be  used 
to  remove  the  second  order  nonlinear  part  of  the  signal.  Since  this  circuit 
uses  diodes,  the  results  of  the  diode  analysis  of  Example  4  can  be  employed 
to  facilitate  this  BDS  analysis.  The  dc-bias  required  to  operate  at  the 
proper  diode  quiescent  point  has  not  been  depicted  since  it  adds  nothing  to 
the  analysis  and  is  blocked  by  the  output  capacitor. 


Fig.  15  Balanced  Diode  Squarer  Circuit 

For  the  purpose  of  analysis,  the  BDS  circuit  can  be  redrawn  as  appears 
in  Fig.  16  The  resistors  R^  and  R^  and  the  resistors  and  have  been 
matched  in  Fig.  16.  Furthermore,  the  two  diodes  are  also  matched;  if  the 
circuit  is  constructed  of  discrete  components,  this  diode  matching  can  be 
accomplished  with  additional  external  resistors  and  capacitors  "hung-on"  the 
diodes. 


The  loops  of  the  circuit  have  been  chosen  so  that  only  a  single  loop 
current  passes  through  a  diode;  this  choice  simplifies  the  algebra  of  the 


analysis.  The  system  or  loop  equations,  and  the  output  equation  for  the 
circuit  of  Fig.  16  can  be  written  as 


1  4 

+Z  I, 

I,. 

a  3 

s  4 

+Z  I, 

I,, 

3  3 

a  4 

y= 

R  Cs+1 

where  Z  =•  — p- —  ,  and  Z  -  Z  +R  +R  .  Now  by  using  the  results  of  the 

S  US  a  S  1  ^ 

diode  analysis  of  Example  2,  the  voltages  u^  and  u^  can  be  written  as  a 
function  of  the  loop  currents.  Then  and  1^^  can  be  expressed  as  Volterra 
NLTFs  operating  on  x.  Then  the  first  three  Volterra  NLTFs  of  the  output  are 
simply  formed. 


Fig.  16  BDS  Circuit  Redrawn 
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where  liberal  operator  notation  is  employed.  There  is  a  similar  expansion 
for  the  second  equation  of  the  above  loop  equations. 

And  finally,  the  powers  of  e  can  be  collected. 
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where  the  diode  Volterra  relationship  is  a  linear  relationship  ;  viz., 
P^(s)=  ^'’(s)=  1/{CqS+X^}. 

In  this  particular  case,  it  is  not  necessary  to  solve  for  the  first 
order  loop  currents  to  obtain  the  first  order  output  y^ .  From  inspection  of 
this  system  of  equations,  because  of  the  symmetry  of  the  equations  it  is 
apparent  that  I  i=  -I  and  I.  =  -I-,,. 


And  therefore. 


Nevertheless,  the  loop  equations  must  be  solved  to  obtain  which  will  be 
required  in  the  evaluation  of  the  higher  order  transfer  functions. 

To  this  end,  equations  (76)  can  be  restated  as 
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It  therefore  follows  that 
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which  can  be  restated  as 
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The  second  order  output  y^  can  now  be  formed  as 
3^2”  ^5^  ^3'^  ^4^ 

But  this  expression  is  unduly  complicated.  To  gain  more  insight  into 

this  second  order  output,  reexamine  Fig.  15.  The  actual  squaring  of  this 

BDS  circuit  is  performed  by  the  resistor-diode  combination  R, "D,  and  f^2~^2' 

The  resistors  and  are  employed  to  sum  the  voltages  u^  and  u^,  and 

therefore  these  two  resistors  usually  have  impedance  values  which  do  not 

"load"  the  preceeding  diode  resistor  combination.  Furthermore,  the 

combination  C  and  R^  are  added  to  eliminate  any  DC-term;  and  it  is  designed 
5 

so  as  not  to  load  the  circuit,  and  to  pass  all  frequencies  of  x.  Therefore, 
the  output  y^  can  be  essentially  described  as  a  filtered  square  term  as 
depicted  in  Fig.  17. 
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^1 

Fig.  17  BDS  Second-Degree  NLTF 


Here  it  is  not  necessary  to  solve  explicitly  for 
Fig.  15;  because  of  the  symmetry  of  the  circuit, 

^1n°  ”^2n 

Now  look  at  the  voltages  u,  and  u„; 


Look  again  at 


^1= 


♦PjCIij] 

*^2^^22^ 

+P3[l23] 

]+P3[-ll3 

]+P^[-I,^]  + 

^P2CIi23 

-P3[I,3] 

So  that  the  voltage  w  and  y  consists  only  of  the  even  orders.  All  odd 
orders  are  zero;  viz., 


y 


3 


=0 


This  same  result  can  be  obtained  by  recognizing  the  symmetry  in  the 
equations  used  to  solve  for  as  was  done  in  solving  for  . 


m 


Summary^. 

As  stated  earlier,  usually  the  effect  of  R_,  R,, ,  R_,  and  C  are  of  no 

i  4  o 

consequences  to  the  application.  But  more  often  the  "filters"  of  and  P 

are  significant  to  the  application.  If  a  cascade  filter  requires  a 
2 

component  KCx  ],  then  the  BDS  circuit  can  be  used  as 
K^[y]-  KCx^] 


where  K  includes  K  operator  and  the  inverse  filter  operations  of  P  and  P_ 
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